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ABSTRACT 

We have compiled a complete extragalactic sample based on ~ 25,000deg^ to a limiting 
flux of 3 X 10~^^ ergscm"^ s"^ 7,000deg2 to a flux limit of 10~" ergs cm"^ s~^) in the 20 - 
40 keV band with INTEGRAL. We have constructed a detailed exposure map to compensate 
for effects of non-uniform exposure. The flux-number relation is best described by a power-law 
with a slope of a = 1.66 ± 0.11. The integration of the cumulative flux per unit area leads to 
/20-40kcV = 2.6 X 10^^" ergs cm^^ s^^ sr^^, which is about 1% of the known 20 — 40keV X-ray 
background. We present the flrst luminosity function of AGN in the 20-40 keV energy range, 
based on 38 extragalactic objects detected by the imager IBIS/ISGRI on-board INTEGRAL. 
The luminosity function shows a smoothly connected two power-law form, with an index of 
7i = 0.8 below, and 72 2.1 above the turn-over luminosity of L» = 2.4 x 10^'^ ergs s^^. The 
emissivity of ah INTEGRAL AGNs per unit volume is W2o-4okcv(> 10''^ ergs s"^) = 2.8 x 
lO^s ergs s~^ h^Q Mpc~^. These results are consistent with those derived in the 2 — 20 keV energy 
band and do not show a signiflcant contribution by Compton-thick objects. Because the sample 
used in this study is truly local {z = 0.022), only limited conclusions can be drawn for the 
evolution of AGNs in this energy band. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies — surveys — galaxies: 
Seyfert 



1. Introduction 

The Galactic X-ray sky is dominated by accret- 
ing binary systems, while the extragalactic sky 
shows mainly active galactic nuclei (AGN) and 
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clusters of galaxies. Studying the population of 
sources in X-ray bands has been a challenge ever 
since the first observations by rocket borne X-ray 
detectors (Giacconi et al. 1962). At soft X-rays 
(0.1-2.4 keV) deep exposures by ROSATYisye re- 
vealed an extragalactic population of mainly broad 
line AGNs, such as type Seyfert 1 and quasars 
(Hasinger et al. 1998; Schmidt et al. 1998). In 
the 2-10 keV range surveys have been carried 
out with ASCA (e.g. Ueda et al. 2001), XMM- 
Newton (e.g. Hasinger 2004), and Chandra (e.g. 
Brandt et al. 2001) and have shown that the 
dominant extragalactic sources are more strongly 
absorbed than those within the ROSAT energy 
band. For a summary on the deep X-ray surveys 
below 10 keV see Brandt & Hasinger (2005). At 



higher energies the data become more scarce. Be- 
tween a few keV and ~ 1 MeV, no all-sky survey 
using imaging instruments has been performed to 
date. The Rossi X-ray Timing Explorer (RXTE) 
sky survey in the 3 — 20 keV energy band revealed 
about 100 AGNs, showing an even higher fraction 
of absorbed {Nh > lO^^cm-^) sources of about 
60% (Sazonov & Revnivtsev 2004). The Interna- 
tional Gamma-Ray Astrophysics Laboratory {IN- 
TEGRAL; Winkler et al. 2003) offers an imprece- 
dented > 20keV collecting area and state-of-the- 
art detector electronics and background rejection 
capabilities. Notably, the imager IBIS with an op- 
erating range from 20 — 1000 keV and a fully-coded 
field of view of 10° x 10° enables us now to study 
a large portion of the sky. A first catalog of AGNs 
showed a similar fraction of absorbed objects as 
the RXTE survey (Beckmann et al. 2006a). The 
Burst Alert Telescope (BAT) of the Swift mission 
(Gehrels et al. 2005) operates in the 15 - 200 keV 
band and uses a detector similar to IBIS/ISGRI, 
but provides a field of view about twice the size. 
The BAT data of the first three months of the 
mission provided a high galactic latitute sur- 
vey, including 44 AGNs (Markwardt et al. 2005). 
Within this sample a weak anti-correlation of lu- 
minosity versus intrinsic absorption was found 
as previously found in the 2 — 10 keV band 
(Ueda et al. 2003; La Franca et al. 2005), reveal- 
ing that most of the objects with luminosities 
Lx > 3 X 10^^ ergs s~^ show no intrinsic absorp- 
tion. Markwardt et al. (2005) also pointed out 
that this luminosity corresponds to the break in 
the luminosity function. 

Related to the compilation of AGN surveys in 
the hard X-rays is the question of what sources 
form the cosmic X-ray background (CXB). While 
the CXB below 20 keV has been the focus of many 
studies, the most reliable measurement in the 10 
- 500 keV has been provided by the High Energy 
Astronomical Observatory {HEAO 1), launched in 
1977 (Marshall et al. 1980). The most precise 
measurement provided by the UCSD/MIT Hard 
X-ray and Gamma- Ray instrument [HEAO 1 A- 
4) shows that the CXB peaks at an energy of 
about 30keV (Marshall et al. 1980, Grubcr et 
al. 1999). The isotropic nature of the X-ray 
background points to an extragalactic origin, and 
as the brightest persistent sources are AGNs, it 
was suggested early on that those objects are the 



main source of the CXB (e.g. Sctti & Woltjer 
1989). In the soft X-rays this concept has been 
proven to be correct through the observations of 
the ROSAT deep X-ray surveys, which showed 
that 90% of the 0.5 - 2.0 keV CXB can be resolved 
into AGNs (Schmidt et al. 1998). At higher ener- 
gies (2 — 10 keV), ASCA and Chandra surveys mea- 
sured the hard X-ray luminosity function (XLF) of 
AGNs and its cosmological evolution. These stud- 
ies show that in this energy range the CXB can 
be explained by AGNs, but with a higher frac- 
tion of absorbed {Nh > lO^^cm"^) objects than 
in the soft X-rays (e.g. Ueda et al. 2003). A 
study based on the RXTE survey by Sazonov & 
Revnivtsev (2004) derived the local hard X-ray lu- 
minosity function of AGNs in the 3 20 kcV band. 
They showed that the summed emissivity of AGNs 
in this energy range is smaller than the total X-ray 
volume emissivity in the local Universe, and sug- 
gested that a comparable X-ray flux may be pro- 
duced together by lower luminosity AGNs, non- 
active galaxies and clusters of galaxies. Using the 
HEAO 1-A2 AGNs, Shinozaki et al. (2006), how- 
ever, obtained a local AGN emissivity which is 
about twice larger than the value of Sazonov & 
Revnivtsev (2004) but consistent with the esti- 
mates by Miyaji et al. (1994) which was based 
on the cross-correlation of the HEAO 1-A2 map 
with IRAS galaxies. 

With the on-going observations of the sky by 
INTEGRAL, a sufficient amount of data is now 
available to derive the AGN hard X-ray lumi- 
nosity function. In this paper we present anal- 
ysis of recent observations performed by the IN- 
TEGRAL satellite, and compare the results with 
previous studies. In Section 2 we describe the 
AGN sample and in Section 3 the methods to de- 
rive the number-flux distribution of INTEGRAL 
AGNs are presented together with the analysis of 
their distribution. Section 4 shows the local lumi- 
nosity function of AGNs as derived from our data, 
followed by a discussion of the results in Section 
5. Throughout this work we applied a cosmology 
with Ho = 70kms-iMpc-i (/170 = 1), fc = 
(flat Universe), flmatter = 0.3, and Aq = 0.7, al- 
though a Ao = and qo = 0.5 cosmology docs not 
change the results significantly because of the low 
redshifts in our sample. 
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2. The INTEGRAL AGN Sample 

Observations in the X-ray to soft gamma-ray 
domain have been performed by the soft gamma- 
ray imager (20-1000 keV) ISGRI (Lebrun et al. 2003) 
on-board the /ATTE'CiJ^L satellite (Winkler et al. 2003). 

The data used here arc taken from orbit revo- 
lutions 19 - 137 and revolutions 142 - 149. The 
list of sources was derived from the analysis as de- 
scribed in Bcckmann ct al. (2006a). The analysis 
was performed using the Offline Science Analy- 
sis (OSA) software version 5.0 distributed by the 
ISDC (Courvoisier ct al. 2003a). Additional ob- 
servations performed later led to further source 
detections within the survey area. We extracted 
spectra at those positions from the data following 
the same procedure. It is understood that most of 
those objects did not result in a significant detec- 
tion > 3(7 in the data set used here, but it ensures 
completeness of the sample at a significance limit 
of 5a (see Section 3). 

The list of 73 sources is shown in Tab. 1. 22 
of the sources have Galactic latitudes — 10° < 6 < 
-1-10° (14, if we only consider the sources with sig- 
nificance > 5cr). In addition to the sample pre- 
sented here, 8 new INTEGRAL sources with no 
identification have been detected in our survey 
with a significance of > Sc. These un-idcntificd 
sources, most of them in the Galactic Plane, are 
not included in this work. The significances listed 
have been derived from the intensity maps pro- 
duced by the OSA software. Different to Beck- 
mann et al. (2006a) we did not use the signifi- 
cances as determined for the whole ISGRI energy 
range by the extraction software, but determined 
the significances based on the count rate and count 
rate error for ISGRI in the 20 40 keV energy 
band only, as this is the relevant energy range 
for this work. Fluxes are determined by integrat- 
ing the best-fit spectral model over the 20 40 kcV 
bandpass. The uncertainty in the absolute flux 
calibration is about 5%. The luminosities listed 
arc the luminosities in this energy band, based 
on the measured (absorbed) flux. The absorp- 
tion listed is the intrinsic absorption in units of 
10^^ cm~^ as measured in soft X-rays below 10 kcV 
by various missions as referenced. We also include 
the most important reference for the INTEGRAL 
data of the particular source in the last coliunn of 
Table 1. The extracted images and source results 



are available in electronic form^. 

In order to provide a complete list of AGNs 
detected by INTEGRAL, we included also those 
sources which are not covered by the data used 
for our study. Those sources are marked in Tab. 1 
and are not used in our analysis. 

3. Number-Flux Distribution of INTE- 
GRAL AGNs 

3.1. Completeness of the Sample 

In order to compute the AGN number-flux re- 
lation it is necessary to have a complete and unbi- 
ased sample. Towards this end, one must under- 
stand the characteristics of the survey, such as the 
sky coverage and completeness for each subset of 
the total sample. Because of the in-homogeneous 
nature of the survey exposure map, we applied a 
significance limit rather than a flux limit to define 
a complete sample. The task is to find a signifi- 
cance limit which ensures that all objects above 
a given flux limit have been included. To test 
for completeness, the Ve/Va-statistic has been ap- 
plied, where Vg stands for the volume that is en- 
closed by the object, and Va is the accessible vol- 
ume, in which the object could have been found 
(Avni & BahcaU 1980). 

In the case of no evolution (Ve/Va) = 0.5 is 
expected. This evolutionary test is applicable only 
to samples complete to a well-defined significance 
limit. It can therefore also be used to test the 
completeness of a sample. We performed a series 
of Ve/Va-tests to the INTEGRAL AGN sample, 
assuming completeness limits in the range of 0. So- 
up to 16(7 ISGRI 20 40 keV significance. For 
a significance limit below the true completeness 
limit of the sample one expects the V^/V^-tests 
to derive a value (Ve/Va) < {Ve/Va)true, where 
{Ve/Va)true is the true test result for a complete 
sample. Above the completeness limit the (Ve/Va) 
values should be distributed around (Ve/Va)true 
within the statistical uncertainties. 

The results of the tests are shown in Figure 1. 
It appears that the sample becomes complete at a 
significance cutoff of approximately 5(7, which in- 
cludes 38 AGNs. The average value is (Vg/Va) = 
0.43 ± 0.05. This is consistent with the expected 
value of 0.5 at the 1.5cr level, suggesting no evo- 

http://heascirc.gsfc.nasa.gov/docs/integral/inthp_Eirchive.html 



3 



lution and a uniform distribution in the local uni- 
verse. It is unlikely that cosmological effects have 
an influence on the result, as the average redshift 
in the sample is z = 0.022, with a maximum red- 
shift of z = 0.13. A positive cosmological evolu- 
tion would result in a slightly higher value than 
0.5. Wc would like to remind that wc use the 
(14/14) test is not to determine any cosmologi- 
cal effects, but use it to see at what significance 
level it returns a stable value. 

3.2. Deriving the Area Corrected Number- 
Flux Distribution 

A correct representation of the number flux dis- 
tribution (i.e. log A^>s versus log 5, see Beckmann 
et al. 2006b) for the sample presented here must 
account for different exposure times comprising 
our survey, and the resulting sensitivity variations. 
We determine here the number density and thus 
the number of AGNs above a given flux has to be 
counted and divided by the sky area in which they 
are detectable throughout the survey. We there- 
fore flrst determined the exposure time in 64, 620 
sky elements of ~ 0.63 deg^ size within our survey. 
In each sky bin, the exposure is the sum of each in- 
dividual exposure multiplied by the fraction of the 
coded field of view in this particular direction. The 
dead time and the good time intervals (GTI) are 
not taken into account but the dead time is fairly 
constant (around 20%) and GTI gaps are very rare 
in IBIS/ISGRI data. Figure 2 shows the exposure 
map in Galactic coordinates for this survey. We 
excluded those fields with an exposure time less 
than 2 ks, resulting in 47,868 sky elements with 
a total coverage of 9.89 sr. The flux limit for a 
given significance limit should be a function of the 
square root of the exposure time, if no systematic 
effects apply, but this assumption cannot be made 
here. The nature of coded-mask imaging leads to 
accumulated systematic effects at longer exposure 
times. In order to achieve a correlation between 
the exposure time and the flux limit, we therefore 
used an empirical approach. For each object we 
computed what we will call its 5(7 equivalent flux 
/scr, based on its actual flux fx and its signiflcance 
s- /sct = fx*5/s. We found a correlation between 
these /5(j values and exposure times, which has a 
scatter of <^ 0.2 dex (Fig. 3). The correlation was 
then fltted by a smooth polynomial of third de- 
gree. This function was then used to estimate the 



limiting flux of each individual survey field. It 
must be noted that the individual limits are not 
important, but only the distribution of those flux 
limits. The total area in the survey for a given 
fiux limit is shown in Figure 4. 

Based on the flux limits for all survey fields, we 
are now able to construct the number flux distri- 
bution for the INTEGRAL AGNs, determining for 
each source flux the total area in which the source 
is detectable with a 5a detection significance in the 
20 — 40 keV energy band. The resulting correlation 
is shown in Figure 5. 

3.3. The Slope of the Number-Flux Dis- 
tribution 

We applied a maximum-likelihood (ML) algo- 
rithm to our empirical number-flux distribution 
to obtain a power-law approximation of the form 
Ni> S) = K ■ S-"'. We note that we are fit- 
ting the "integrated" A^(> S) function, as distinct 
from the " differential" number-flux function. The 
latter entails binning the data, and thus some loss 
of information is incurred. The advantage of flt- 
ting the differential distribution is that a simple 
least squares procedure may be employed. How- 
ever, given the modest size of our sample, the ex- 
pected loss of accuracy was considered unaccept- 
able. 

Our approach was based on the formalism de- 
rived by Murdoch, Crawford & Jauncey (1973), 
also following the implementation of Piccinotti 
et al. (1982). The latter involved modiflcation 
of the basic ML incorporated to facilitate han- 
dling of individual source flux-measurement er- 
rors. The ML method also involves the applica- 
tion of a Kolmogorov-Smirnov (K-S) test as part 
of the procedure to optimize the fit, as detailed 
in Murdoch, Crawford & Jauncey (1973) (we note 
that the KS test as applied in this context is not a 
measure of the overall goodness of fit). Once the 
slope is determined, a chi-square minimization is 
used to determine the amplitude K. 

For this analysis, we used the complete sub- 
sample of 38 sources for which the statisti- 
cal significance of our flux determinations was 
at a level of 5a or greater. The dimmest 
source among this sub-sample was fx = 5.6 x 
10""'^^ ergs cm^^ s^^, and the brightest was fx = 
3.2 X 10"^°ergscm-2s"^ We derived a ML 
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probability distribution, which can be approxi- 
mated by a Gaussian, with our best fit param- 
eters of a = 1.66 ±0.11. A normalization of 
K = 0.44 sr"^ (10"^" ergs cm^^s^i)" was then 
obtained by performing a least-squares fit, with 
the slope fixed to the ML value. This calculation 
did not take into account possible inaccuracies as- 
sociated with scatter in Fig. 3, and thus in the 
true detection limit. The true exposure time is 
also affected by small variations in dead-time ef- 
fects known to occur in the ISGRI detector. A 
conservative upper limit on the exposure time un- 
certainty is < 2%. This leads to uncertainty in 
the final log N - log S primarily manifest in the 
normalization and it should not affect the slope 
significantly. Furthermore, the uncertainty in the 
detection limit will affect mainly the low flux end 
of the Log(N) - Log(S) distribution. The high flux 
end is less sensitive to scatter, since it is based 
on a larger sky area (Fig. 4). To make a more 
quantitative assessment, we have recomputed the 
ML Log(N) - Log(S) calculation for scenarios in 
which the exposure time - flux limit curve shifted 
in amplitude and pivoted about the 700 ks point 
where we have the highest density of measure- 
ments. For those scenarios, we found that the 
inferred Log(N) - Log(S) slope varied by less than 
about 5%, which is contained within the range of 
our quoted 1-sigma uncertainty. The amplitude 
varied by as much as 7% in the extreme case, but 
for the pivoted cases, by only a few percent. We 
thus conclude that the maximum uncertainty re- 
sulting from possible systematics in our effective 
area correction is bounded by about 5% in slope 
and 7% in amplitude. 

4. The Local Luminosity Function of AGNs 
at 20 ~ 40 keV 

The complete sample of INTEGRAL AGNs 
with a detection significance > 5a also allows us 
to derive the density of these objects in the local 
Universe as a function of their luminosity. In or- 
der to derive the density of objects above a given 
luminosity, one has to determine for each source in 
a complete sample the space volume in which this 
source could have been found considering both the 
flux limit of each survey field and the flux of the 
object. We have again used the correlation be- 
tween exposure time and flux limit as discussed 
in the previous section in order to assign a 5a flux 



limit to each survey field. Then the maximum red- 
shift 

Zmax s-t which an object with luminosity Lx 
would have been detectable in each sky element 
was used to compute the total accessible volume 



^ n 



(1) 



with A'' being the number of sky elements in which 
the object would have been detectable and fli the 
solid angle covered by sky element i, and Vi the 
enclosed volume based on the maximum redshift 
at which the object could have been detected in 
this sky element. Figure 6 shows the cumula- 
tive luminosity function for 38 INTEGRAL de- 
tected (> 5a) AGNs in the 20 - 40 keV energy 
band. Here the density (j) describes the number 
of objects per Mpc^ above a given luminosity Lx- 

K 

(j) = ^ V~l with K being the number of objects 
1=1 

with luminosities > Lx- Blazars have been ex- 
cluded because their emission is not isotropic. The 
redshifts in the sample range from z = 0.001 to 
z = 0.129 with an average redshift of ^ = 0.022. 
Thus the luminosity function is truly a local one. 
Figure 7 shows the luminosity function in differ- 
ential form. In this presentation the data points 
arc independent of each other. In case one of the 
luminosity bins would suffer from incompleteness 
compared to the other bins, this would result in a 
break or dip in the differential luminosity function. 
The errors are based on the number of objects con- 
tributing to each value. The differential XLF also 
shows, like the cumulative one, a turnover around 
Lx = {5- 10) X 10^3 ergs s-i. 

Because our study is based solely on low red- 
shift objects, we are not able to constrain models 
involving evolution with redshift. Nevertheless we 
can compare the XLF presented here with model 
predictions from previous investigations. XLFs 
are often fit by a smoothly connected two power- 
law function of the form (Maccacaro et al. 1991) 



d(j>{Lx,z = 0) 



dlogL 



X 



h2L 

Ldf. 



71 



+ 



72 



Lx 

(2) 

We fit this function using a least-squares 
method applying the Levenberg-Marquardt al- 
gorithm (Marquardt 1963). The best fit values 
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wc obtained arc A ^ x 10~^ h^^Mpc'^ , 

71 = 0.80 ± 0.15, 72 = 2.11 ± 0.22, and logL* = 
43.38 ± 0.35 with in units of h^Q eTgss~^. 
The la errors have been determined by applying 
a Monte Carlo simulation which simultaneously 
takes into account the flux errors on the indi- 
vidual sources, the error induced by deriving an 
average luminosity per bin, and the statistical er- 
ror of the density based on the number of objects 
contributing to the density value. Each simulated 
data set included 9 luminosity values with a den- 
sity value for each of them. These values where 
then fit by the smoothly connected two power-law 
function as described above. The scatter in the 
resulting parameters gave the error estimates as 
shown above. 

The parameter values describing the differen- 
tial luminosity function are consistent with val- 
ues derived from the 2-10 keV XLF of AGNs 
as shown by e.g. Ueda et al. (2003), Franca 
et al. (2005), and Shinozaki et al. (2006). For 
example the work by Ueda et al. (2003) reveals 
for a pure density evolution model the same val- 
ues (within the error bars) for 71 and 72, but a 
higher logL, = 44.11 ± 0.23. The higher value 
can be easily explained by the different energy 
bands applied. A single power law with photon 
index of F = 2 in the range 2 — 40keV would 
lead to i(2_iokeV)/-i'(20-40keV) = 2.3, assuming 
no intrinsic absorption. This has, of course, no 
implications for the XLF at higher redshifts. The 
values are also consistent with the luminosity func- 
tion for AGNs in the 3 — 20 keV band as derived 
by Sazonov & Revnivtsev (2004) from the RXTE 
all-sky survey. 

Information about intrinsic absorption is avail- 
able for 32 of the 38 objects (89%) from soft X-ray 
observations. This enables us to derive the lumi- 
nosity function for absorbed (Nh > 10^^ cm^^) 
and unabsorbed sources, as shown in Figure 8. 
The absorbed sources have a higher density than 
the unabsorbed sources at low luminosities, while 
this trend is inverted at high luminosities. The 
luminosity where both AGN types have similar 
densities is about i(2o-40kcV) = 3 x 10"''' ergs^^. 
This tendency is also evident when comparing the 
fraction of absorbed AGNs with the luminosity 
in the three luminosity bins depicted in Figure 9. 
The luminosity bins have been chosen so that an 
equal number of objects are contained in each bin. 



The position of the data point along the luminos- 
ity axis indicates the average luminosity in this 
bin, while the error bars in luminosity indicate 
the range of luminosities covered. A comparable 
trend has been seen also below 10 keV. In a re- 
cent study of HEAO-1 data of 49 AGNs, Shinozaki 
et al. (2006) showed that the XLF for absorbed 
AGNs drops more rapidly (72 = 3.34^°;^°) at 
higher luminosities than that of unabsorbed AGNs 
(72=2.34;°-^). 

Based on the luminosity function, the contribu- 
tion of the AGNs to the total X-ray emissivity W 
can be estimated (Sazonov & Revnivtsev 2004). 
This can be done by simply multiplying the XLF 
by the luminosity in each bin and integrating 
over the range of luminosities (lO"'^ ergs s"'^ < 
-^20-40 keV < lO*^'^ crgs s~^) . This results in 
W2o-40kcv(> 10^1 ergs s-1) = (2.8 ± 0.8) x 
10"^* ergs s^^ h^Q Mpc""^. Please note that absorp- 
tion does not affect the luminosities in this energy 
range and therefore the values given here are in- 
trinsic emissivities. 

5. Discussion 

A simple power-law model fitted to the num- 
ber flux distribution (Fig. 5) has a slope of a = 
1.66 ± 0.11. Even though the difference from the 
Eucledian value is not statistically significant, at 
a 1.5cr level, a deviation from this value could 
have two reasons. The difference might indicate 
that the area density at the low flux end of the 
distribution has been slightly over corrected. One 
has to keep in mind that only a few sources de- 
rived from a small area of the sky are constrain- 
ing the low flux end. Another reason for the dif- 
ference could be that the distribution of AGNs 
in the very local universe is not isotropic, caused 
e.g. by the local group and other clustering of 
galaxies. Krivonos et al. (2005) studied the ex- 
tragalactic source counts as observed by INTE- 
GRAL in the 20-50 keV energy band in the Coma 
region. Based on 12 source detections they deter- 
mine a surface density of (1.4 ± 0.5) x 10~^ deg~^ 
above a threshold of 10^^^ ergscm"^ s~^ in the 
20 — 50 keV energy band, where we get a consis- 
tent value of (1.2 ±0.2) x 10~^deg~^. Comparing 
the total flux of all the objects in the AGN sample 
(.f20-40kcv = 2.6 X 10^-^° ergscm"^ s^^ sr~^) with 
the flux of the X-ray background as presented by 
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Gruber et al. (1999) shows that the INTEGRAL 
AGN account only for about 1% of the expected 
value. This is expected when taking into account 
the high flux limit of our sample: La Franca et al. 
(2005) have shown that objects with /2-iokeV > 
10~^^ ergs cm~^ contribute less than 1% to the 
CXB in the 2 — lOkeV energy range. This flux 
limit extrapolates to the faintest flux in our sam- 
ple of /2o-40keV = 5.6 X 10~^^ ergs cm~^ s~^ for a 
r = 1.9 power law spectrum. 

We compared the unabsorbed emissivity per 
unit volume of our objects W2o-40kov(> 10'*^ ergs s^^) = 
2.8 X 10^® ergss"^ hfo Mpc~^ with that observed 
by RXTE in the 3-20 keV band. Assuming 
an average power law of F = 2, the extrap- 
olated value is W3_20kev(> 10^^ ergss"-^) = 
(7.7 ± 2.2) X lO^Sergss-^hfoMpc-^, which is a 
factor of 2 larger than the value measured by 
RXTE (Sazonov & Revnivtsev 2004) but consis- 
tent within the la error. If we apply the con- 
version to the 2 — lOkeV energy band, we de- 
rive the intrinsic emissivity W2-iokcV = (6.4 ± 
1.8) X 10^^ergss~^Mpc~^, consistent with the 
value derived from the HEAO-1 measurements 
{W = (5.9 ±1.2) X 10^8 ergs s-iMpe-3; Shinozaki 
et al. 2006). This showed that the local X-ray 
volume emissivity in the 2-10 keV band is consis- 
tent with the emissivity from AGNs alone. It has 
to be pointed out that the value derived from our 
sample and the one based on HEAO-1 data are 
higher than the one based on the RXTE All-Sky 
Survey {W = (2.7 ± 0.7) x lO^^ ergss"^ Mpc'^; 
Sazonov & Revnivtsev 2004). 

The luminosity function derived from the IN- 
TEGRAL 20 - 40 keV AGN sample appears to be 
consistent with the XLF in the 2 — 20 keV range. 
A turnover in the XLF at ~ 2.4 x 10''^ ergs s"^ 
is observed (Fig. 7). Below this luminosity also 
the fraction of absorbed AGNs starts to exceed 
that of the unabsorbed ones, although the ef- 
fect is significant at only on a la level (Fig. 8). 
Both efl:ects have been seen also in the 2 — 
10 keV (Ueda et al. 2003; La Franca et al. 2005; 
Shinozaki et al. 2006) and in the 3 — 20 keV band 
(Sazonov & Revnivtsev 2004). This implies that 
we do detect a similar source population as at 
lower energies. 

If a larger fraction of absorbed AGNs is neces- 
sary to explain the cosmic X-ray background at ~ 
30 keV as indicated by HE AO 1 A-4 measurements 



(Gruber et al. 1999), the fraction of absorbed 
sources could be correlated with redshift. It has 
for example been proposed that there is an evolu- 
tion of the population leading to a higher fraction 
of absorbed sources at higher redshifts. It should 
be noted however that this effect is not clearly 
detectable in the 2 — 10 keV range. The fraction 
of absorbed sources seems to depend on luminos- 
ity (Ueda et al. 2003; Treister & Urry 2005), as is 
also seen in the 20 - 40keV band (Fig. 9). But 
some studies come to the conclusion that there 
is no evolution of intrinsic Nh (Ueda et al. 2003; 
Treister & Urry 2005), while others find the frac- 
tion of absorbed sources increasing with redshift 
(La Franca et al. 2005). La Franca et al. also 
find that a combination of effects (the fraction of 
absorbed AGNs decreases with the intrinsic X-ray 
luminosity, and increases with the redshift) can 
be explained by a luminosity-dependent density 
evolution model. They further show that the lu- 
minosity function of AGNs with low luminosities 
as those presented here peaks at z ^ 0.7 while high 
luminosity AGNs peak at z ~ 2. Unified models 
also predict, depending on the applied model, a 
fraction of absorbed AGNs of 0.6 — 0.7 compared 
to the total population for high-flux low-redshift 
objects (Treister & Urry 2005). Worsley et al. 
(2005) examined Ghandra and XMM-Newton deep 
fields and come to the conclusion that the miss- 
ing CXB component is formed by highly obscured 
AGNs at redshifts ^ 0.5 — 1.5 with column densi- 
ties of the order of fx = lO^^ - lO^" cm'^. Evi- 
dence for this scenario is also found in a study of 
Chandra and Spitzer data (PoUetta et al. 2006). 
Combining multiwavelength data, this work esti- 
mates a surface density of 25AGNdeg~^ in the 
infrared in the 0.6 deg^ Chandra /SWIRE field, 
and only 33% of them are detected in the X-rays 
down to /o.3-8keV = 10~^^ crgs cm~^ s~^. The 
work also indicates a higher abundance of lumi- 
nous and Compton-thick AGNs at higher red- 
shifts {z » 0.5). This source population would 
be missed by the study presented here, because of 
the low redshifts {z = 0.022) of the INTEGRAL 
AGNs. 

Several studies (Ueda et al. 2003; Treister & Urry 2005) 
propose that the absorbed AGNs needed to ex- 
plain the CXB should be Compton thick, and 
therefore would have been missed at 2 — lOkeV. 
This argument does not hold for the INTEGRAL 
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observations, where the impact of absorption is 
much less severe than at lower energies. The ef- 
fect on the measured flux of a source with pho- 
ton index r = 2 for Compton thick absorption 
{Nh = 10^'* cm~^) is only a 5% decrease in flux 
(40% for Nh = lO^^cm-^). It is therefore un- 
likely that many Compton-thick objects have been 
missed by the INTEGRAL studies performed to 
date. One possibility would be, that they are 
among the newly discovered sources found by IN- 
TEGRAL. The fraction of unidentifled objects 
among the INTEGRAL discovered sources is ap- 
proximately 50%. Eight such sources without 
cross-identiflcation have a signiflcance above Sa- 
in the data set discussed here. Thus, if they are 
ultimately identified as AGNs, they would have to 
be considered in this study. It should be pointed 
out though, that most of the sources discovered 
by INTEGRAL are located close to the Galactic 
plane and are more likely to belong to the Galaxy: 
the Second IBIS/ISGRI Soft Gamma-Ray Survey 
Catalog (Bird et al. 2006) lists 55 new sources 
detected by INTEGRAL, of which 93% are lo- 
cated within -10° < 6 < -|-10°. Among these 
55 sources, 3 are listed as extragalactic sources, 
18 are of Galactic origin, and 29 have not been 
identified yet. 

In addition, those objects which have been clas- 
sified as AGNs based on soft X-ray and/or op- 
tical follow-up studies, are no more likely to be 
Compton-thick objects than the overall AGN pop- 
ulation studied here. Only four AGNs (NGC 1068, 
NGC 4945, MRK 3, Circinus galaxy) detected 
by INTEGRAL have; been proven to be Comp- 
ton thick objects so far, and none of them showed 
absorbtion of Nh > 5 x lO^^cm"^. In order to 
clarify this point, observations at soft X-rays of 
those objects without information about intrinsic 
absorption are required for all INTEGRAL de- 
tected AGNs (Tab. 1). At present 23 % of the 
INTEGRAL AGN are missing absorption infor- 
mation. A first indication of what the absorp- 
tion in these sources might be, can be derived 
from comparison of the INTEGRAL fluxes with 
ROSAT All-Sky Survey (RASS) Faint Source Cat- 
alogue data (Voges et al. 2000). In order to do so 
we assumed a simple power law with photon in- 
dex r = 2.0 between the ROSAT 0.1 - 2.4 keV 
band and the INTEGRAL 20 - 40 keV range and 
flt the absorption. In the six cases where no de- 



tection was achieved in the RASS, an upper limit 
of /(o.i-2.4keV) < 10"^'^ crgs cm~^ s~^ has been as- 
sumed, resulting in a lower limit for the absorption 

Nh > {5- 11) X 10^2 cm"2 pjg iq gj^g^ 

the distribution of intrinsic absorption. It has to 
be pointed out that the estimated values can only 
give an idea about the distribution of intrinsic ab- 
sorption and should not be taken literally, as the 
spectral slope between the measurements is un- 
known and the observations are not simultane- 
ous. Nevertheless apparently none of the RASS 
detections and non-detections requires an intrin- 
sic absorption of Nh > 2 x 10^"^cm~^. Therefore 
it appears unlikely that a significant fraction of 
INTEGRAL AGNs wifl show an intrinsic absorp- 
tion Nh > lO^'^cm^^. However, if we assume 
that the RASS non-detections are all Compton 
thick AGNs, the fraction of this class of sources 
rises from 6% to 16% when considering all 63 non- 
blazar AGNs seen by INTEGRAL, and from 8% to 
13% for the complete sample with 38 objects. This 
is in good agreement with the fraction of 11% of 
Compton thick AGN as seen in the Swift/BAT sur- 
vey (Markwardt et al. 2005). The picture is less 
clear when referring to the optical classification. 
Here the INTEGRAL survey finds 12 Seyfert 1 
(33%), 14 Seyfert 2, and 10 intermediate Seyfert 
1.5 in the complete sample, while the Swift/BAT 
survey contains only 20% of type 1 Seyfert galax- 
ies. It should be pointed out though that the 
classification based on intrinsic absorption gives 
a more objective criterion in order to define AGN 
subclasses than the optical classification with its 
many subtypes. The finding of the BAT survey 
that virtually all sources with logLx < 43.5 are 
absorbed, cannot be confirmed by our study, in 
which we detect a fraction of 33% of sources with 
Nh < 10^^ cm~^ among the sources below this 
luminosity. This also reflected in the observation 
that although the absorbed sources become more 
dominant below this luminosity, the trend is not 
overwhelmingly strong (Fig. 8). 

Most investigations to date have been focused 
on the X-rays below 20 keV, and INTEGRAL 
can add substantial information to the nature of 
bright AGNs in the local Universe. Consider- 
ing the expected composition of the hard X-ray 
background, it does not currently appear that 
the population detected by INTEGRAL can ex- 
plain the peak at 30 keV, as Compton thick 
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AGNs arc apparently less abundant than ex- 
pected (Treister & Urry 2005). But this picture 
might change if we assume that all INTEGRAL 
AGNs lacking soft X-ray data and without counter 
parts in the RASS to be Compton thick. In ad- 
dition the sample presented here might be still 
too small to constrain the fraction of obscured 
sources, and the missing Compton thick AGNs 
could be detectable when studying sources with 
/(20-40keV) < 10"^^ergscm-2s-^ 

6. Conclusions 

The extragalactic sample derived from the IN- 
TEGRAL public data archive comprises 63 low 
redshift Seyfert galaxies {{z) = 0.022 ±0.003) and 
8 blazars in the hard X-ray domain. Two galaxy 
clusters are also detected, but no star-burst galaxy 
has been as yet. This INTEGRAL AGN sample 
is thus the largest one presented so far. 38 of the 
Seyfert galaxies form a complete sample with sig- 
nificance limit of 5(7. 

The number flux distribution is approximated 
by a power-law with a slope of a = 1.66 ± 0.11. 
Because of the high flux limit of our sample the 
objects account in total for less than 1% of the 
20 — 40 keV cosmic X-ray background. The emis- 
sivity of all AGNs per unit volume Vr20-40kov(> 
lO-^^ergss-i) = 2.8 x lO^^ ergss'^ h^o Mpc'^ ap- 
pears to be consistent with the background es- 
timates in the 2-10 keV energy band based on 
the cross-correlation of the HE AO 1-A2 map with 
/R^ 5 galaxies (Miyaji et al. 1994). 

The luminosity function in the 20 — 40 kcV en- 
ergy range is consistent with that measured in the 
2 — 20keV band. Below the turnover luminosity 
of = 2.4 X 10'*-''crgss-i the absorbed AGNs 
become dominant over the unabsorbed ones. The 
fraction of Compton thick AGNs with known in- 
trinsic absorption is found to be small (8%) in 
our AGN sample. For the sources without reli- 
able absorption information we derived an esti- 
mate from the comparison with ROSAT All-Sky 
Survey data and find that the data do not re- 
quire additional Compton thick objects within the 
sample presented here. It has to be pointed out 
though, that the sources without RASS counter- 
part could be Compton thick which would increase 
the ratio of this source type to 13% in the com- 
plete sample. Evolution of the source population 



can play a major role in the sense that the frac- 
tion of absorbed sources among AGNs might be 
correlated with redshift, as proposed for example 
by Worslcy et al. (2005). 

Over the life time of the INTEGRAL mission 
we expect to detect of the order of 200 AGNs. 
Combining these data with the studies based on 
Swift/BAT, operating in a similar energy band as 
IBIS/ISGRI, will further constrain the hard X- 
ray luminosity function of AGNs. But we will 
still be limited to the relatively high fltix end of 
the distribution. Because of this INTEGRAL and 
Swift/BAT will most likely not be able to test evo- 
lutionary scenarios of AGNs and thus will be inad- 
equate to explain the cosmic X-ray background at 
E > 20keV. Future missions with larger collect- 
ing areas and/or focusing optics will be required 
to answer the question of what dominates the Uni- 
verse in the hard X-rays. 
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Fig. 1. — Result of the Ve/Va test, assuming differ- 
ent completeness limits in ISGRI significance. The 
dot-dashed line shows the average (Ve/Va) value 
for the objects with significance > 5a. 




Fig. 2. — Exposure map representing the data 

used in our analysis in Galactic coordinates. Con- 
tours indicate 2 ks, 5 ks, 10 ks, 20 ks, 100 ks, 
200 ks, 500 ks, 1 Ms, 2 Ms exposure time. INTE- 
GRAL spent most of the observing time on and 
near the Galactic plane, with a strong focus on the 
Galactic center and on areas including bright hard 
X-ray sources like the Cygnus region. Vela, GRS 
1915+105, and the Crab. Fields at high Galac- 
tic latitude include Coma, Virgo and the region 
around NGC 4151. 
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Table 1 
INTEGBAL AGN cataloc 



Name 


Type 


z 


R.A. 




Dccl. 




a 

cxp. 


ISGRI 


f b 

JX 


log L20-40 kcV 




Kei. 








(J2000.0) 


(J2000.0) 


(ks) 






(ergs s~i) 






lES 0033+595 


BL Lac 


0.086 


00 


35 


53 


+59 


50 


05 


1449 


3.5 


0.37 


43.83 


0.36^ 


2 


NGC 788 


Sy 1/2 


0.0136 


02 


01 


06 


-06 


48 


56 


311 


10.7 


2.98 


43.09 


< 0.02^ 


4 


IGR J02097+5222 


Sy 1 


0.0492 


02 


09 


46 


+52 


22 


48 


26 


4.9 


3.8 


44.34 


? 


30 


NGC 1068 


Sy 2 


0.003793 


02 


42 


41 


-00 


00 


48 


311 


4.3 


0.93 


41.47 


> 150^ 


4 


QSO B0241+62 


Sy 1 


0.044557 


02 


44 


58 


+62 


28 


07 


43 


3.4 


2.02 


43.97 


1.58 


7 


NGC 1142 


Sy 2 


0.028847 


02 


55 


12 


-00 


11 


01 


311 


5.5 


1.58 


43.48 


? 


8 


NGC 1275 


Sy 2 


0.017559 


03 


19 


48 


+41 


30 


42 


264 


8.4 


1.89 


43.12 


3.75'' 


4 


3C 111 


Sy 1 


0.048500 


04 


18 


21 


+38 


01 


36 


67 


10.0 


6.27 


44.54 


0.63^ 


4 


UGC 3142 


Sy 1 


0.021655 


04 


43 


47 


+28 


58 


19 


247 


16.8 


5.46 


43.76 


7 


2 


LEDA 168563 


Sy 1 


0.0290 


04 


52 


05 


+49 


32 


45 


28 


2.8 


2.27 


43.64 


7 


2 


MCG +8-11-11 


Sy 1.5 


0.020484 


05 


54 


54 


+46 


26 


22 


21 


6.2 


6.07 


43.76 


< 0.02*' 


4 


MRK 3 


Sy 2 


0.013509 


06 


15 


36 


+71 


02 


15 


472 


15.9 


3.65 


43.17 


110" 


4 


MRK 6 


Sy 1.5 


0.018813 


06 


52 


12 


+74 


25 


37 


482 


8.7 


2.01 


43.21 


lO** 


4 


S5 0716+714 


BL Lac 


0.3'^ 


07 


21 


53 


+71 


20 


36 


482 


0.7 


0.14 


44.41'^ 


< 0.01=* 


9 


ESQ 209-12 


Sy 1.5 


0.040495 


08 


01 


58 


-49 


46 


36 


1543 


6.7 


0.86 


43.52 


? 


7 


FRL 1146 


Sy 1 


0.031789 


08 


38 


31 


-35 


59 


35 


849 


3.6 


0.60 


43.15 


7 


7 


S5 0836+710 


FSRQ 


2.172 


08 


41 


24 


+70 


53 


42 


391 


6.4 


1.73 


47.79 


0.11^ 


9 


MCG-05-23-16 


Syl.9 


0.008486 


09 


47 


40 


-30 


56 


56 


2 


2.3 


11.20 


43.25 


1.6i« 


16 


IGR J10404-4625 


Sy 2 


0.0237 


10 


40 


22 


-46 


25 


26 


46 


1.5 


0.67 


42.93 


7 


10 


NGC 3783 


Sy 1 


0.00973 


11 


39 


02 


-37 


44 


19 


18° 


5.6 


6.2 


43.11 


< 0.4^ 


23 


IGR J12026-5349 


AGN 


0.028 


12 


02 


48 


-53 


50 


08 


191 


5.5 


1.86 


43.52 


2.2^1 


11 


NGC 4051 


Sy 1.5 


0.002336 


12 


03 


10 


+44 


31 


53 


443 


8.4 


1.80 


41.34 


< O.Ol'^ 


4 


NGC 4151 


Sy 1.5 


0.003320 


12 


10 


33 


+39 


24 


21 


483 


163.3 


26.13 


42.80 


6.9^^ 


12 


NGC 4253 


Sy 1.5 


0.012929 


12 


18 


27 


+29 


48 


46 


715 


6.1 


0.93 


42.54 


0.8*' 


4 


4C +04.42 


BL Lac 


0.965 


12 


22 


23 


+04 


13 


16 


690 


4.5 


0.80 


46.58 


7 


7 


NGC 4388 


Sy 2 


0.008419 


12 


25 


47 


+ 12 


39 


44 


215 


34.8 


9.54 


43.18 


27I3 


13 


NGC 4395 


Sy 1.8 


0.001064 


12 


25 


49 


+33 


32 


48 


739 


5.1 


0.56 


40.14 


0.15^ 


4 


3C 273 


Blazar 


0.15834 


12 


29 


07 


+02 


03 


09 


655 


34.2 


5.50 


45.58 


0.5* 


14 


NGC 4507 


Sy 2 


0.011801 


12 


35 


37 


-39 


54 


33 


152 


14.9 


6.46 


43.30 


29** 


4 


IGR J12391-1612 


Sy 2 


0.036 


12 


39 


06 


-16 


10 


47 


83 


1.4 


3.46 


44.02 


1.9" 


11 


NGC 4593 


Sy 1 


0.009000 


12 


39 


39 


-05 


20 


39 


723 


20.1 


3.31 


42.78 


0.02** 


4 


IGR J12415-5750 


Sy 2 


0.024 


12 


41 


24 


-57 


50 


24 


440 


1.1 


0.33 


42.64 


7 


15 


3C 279 


Blazar 


0.53620 


12 


56 


11 


-05 


47 


22 


497 


3.6 


0.82 


45.97 


< 0.13^ 


4 


Coma Cluster 


GClstr 


0.023100 


12 


59 


48 


+27 


58 


48 


516 


7.2 


1.09 


43.11 


< 0.01* 


4 


NGC 4945 


Sy 2 


0.001878 


13 


05 


27 


-49 


28 


06 


276 


33.8 


9.85 


41.88 


400" 


16 


ESQ 323-G077 


Sy 2 


0.015014 


13 


06 


26 


-40 


24 


53 


761 


6.9 


1.20 


42.78 


551^ 


17 


IGR J13091+1137 


AGN 


0.0251 


13 


09 


06 


+ 11 


38 


03 


48 


2.0 


1.06 


43.18 


90" 


15 


NGC 5033 


Sy 1.9 


0.002919 


13 


13 


28 


+36 


35 


38 


377 


4.6 


1.06 


41.30 


2.9'' 


4 


Cen A 


Sy 2 


0.001830 


13 


25 


28 


-43 


01 


09 


532 


167.4 


32.28 


42.38 


12.5* 


16 


MCG-06-30-015 


Sy 1 


0.007749 


13 


35 


54 


-34 


17 


43 


567 


4.9 


0.73 


41.99 


7.7'^ 


4 


4U 1344-60 


Sy 1.5 


0.012 


13 


47 


25 


-60 


38 


36 


603 


16.6 


2.83 


43.02 


5^3 


4 


IC 4329A 


Sy 1.2 


0.016054 


13 


49 


19 


-30 


18 


36 


440 


41.7 


8.19 


43.68 


0.42'' 


4 


Circinus gal. 


Sy 2 


0.001448 


14 


13 


10 


-65 


20 


21 


589 


58.9 


10.73 


41.69 


360'' 


16 


NGC 5506 


Sy 1.9 


0.006181 


14 


13 


15 


-03 


12 


27 


101 


12.6 


4.21 


42.55 


3.415 


4 


NGC 5548 


Sy 1.5 


0.017175 


14 


18 


00 


+25 


08 


12 


211° 


2.4 


0.71 


42.67 


0.51'^ 


4 


PG 1416-129 


Sy 1 


0.129280 


14 


19 


04 


-13 


10 


44 


117 


8.3 


4.86 


45.33 


0.09* 


4 


ESQ 511-G030 


Sy 1 


0.022389 


14 


19 


22 


-26 


38 


41 


145 


5.1 


1.93 


43.34 


< 0.05'^ 


17 


IC 4518 


Sy 2 


0.015728 


14 


57 


43 


-43 


07 


54 


338 


2.3 


0.49 


42.44 


7 


4 


IGR J16119-6036 


Sy 1 


0.016 


16 


11 


54 


-60 


36 


00 


475 


1.2 


0.25 


42.16 


7 


2 


IGR J16482-3036 


Sy 1 


0.0313 


16 


48 


17 


-30 


35 


08 


973 


4.2 


0.73 


43.22 


7 


10 


NGC 6221 


Sy 1/2 


0.004977 


16 


52 


46 


-59 


13 


07 


523 


5.6 


1.32 


41.86 


^18 


4 


Oph Cluster 


GClstr 


0.028 


17 


12 


26 


-23 


22 


33 


1763 


30.8 


4.10 


43.90 


7 


25 


NGC 6300 


Sy 2 


0.003699 


17 


17 


00 


-62 


49 


14 


173 


10.0 


3.91 


42.07 


22I3 


16 


GRS 1734-292 


Sy 1 


0.021400 


17 


37 


24 


-29 


10 


48 


3332 


45.9 


4.03 


43.62 


3.7* 


24 


2E 1739.1-1210 


Sy 1 


0.037 


17 


41 


54 


-12 


11 


52 


631 


5.5 


1.03 


43.51 


7 


7 


IGR J18027-1455 


Sy 1 


0.035000 


18 


02 


47 


-14 


54 


55 


942 


12.6 


2.03 


43.76 


19.0* 


20 


PKS 1830-211 


Blazar 


2.507 


18 


33 


40 


-21 


03 


40 


1069 


12.7 


2.07 


48.02 


< 0.7^ 


21 



12 




-11,5 

4.5 5 5.5 6 

log t. [sec] 

Fig. 3. — Correlation of exposure time and flux 
limit (5cr, 20 - 40keV) for the AGNs in this study 
(see Sect. 3. 3). The curve shows a smooth poly- 
nomial fit for flux limit versus logarithm of the 
exposure time. 
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log fj.^.^ [ergs cm"^ s"^] 

Fig. 4. — The total survey area with respect to 
the 5a flux limit in the 20 - 40 keV. The curve 
is based on all 47,868 sky elements of the survey 
with an exposure of at least 2 ks. For a flux limit 
of /jsf > 3 X 10~^^ ergscm"^ the survey covers 
76% of the sky, and for fx > 10~^^ ergscm~^s~^ 
17%. 



Fig. 5. — Number flux distribution per steradian 
of INTEGRAL AGNs with a detection significance 
> 5a. Blazars have been excluded. The maxi- 
mum likelihood slope as described in Section 3.3 
is 1.66 ±0.11. 
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Fig. 6. — Cumulative luminosity function of IN- 
TEGRAL AGNs with a detection significance > 
5a. Blazars have been excluded. The density (f> 
describes the number of objects per Mpc^ above 
a given luminosity Lx. As an example some er- 
ror bars indicating the Poissonian error are shown. 
Note that the errors are not independent of each 
other. 
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Fig. 7. — Differential luminosity function of AGNs 
with AlogLx = 0.5 binning. The line shows a 
fit to a smoothly connected two power-law func- 
tion with a turnover luminosity at — 2.4 x 
10"'^ ergss^^ 




41 42 43 44 45 

log 1-20-40 kev [ergs/sec] 

Fig. 8. — Cumulative AGN luminosity function 
for 19 absorbed {Nh > lO^^cm"^; triangles) and 
12 unabsorbed sources (octagons). As an example 
some error bars indicating the Poissonian error are 
shown. 
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Fig. 9.~ Fraction of absorbed AGNs {Nh > 
10^^ cm~^) versus luminosity. The position of the 
data point along the luminosity axis indicates the 
average luminosity in this bin, while the error bars 
in luminosity indicate the range of luminosities 
covered. 



14 



10 - 




80 22 84 



log N|, I cm ^ I 

Fig. 10. — Distribution of intrinsic absorption for 
all INTEGRAL AGNs (blazars excluded), as mea- 
sured in the soft X-rays. The shaded area shows 
the reliable measurements, the other values are 
based on comparison of ROSAT All-Sky Survey 
and INTEGRAL data. Lower limits on absorp- 
tion have been excluded. 
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Table 1 — Continued 



Name Type z R.A. Decl. exp." ISGRI fx^ logi20-40kev Nh°'^ Ref/ 

(J2000.0) (J2000.0) (ks) (cr) (ergs s-^) 



ESQ 103-G35 


Sy 2 


0, 


.013286 


18 


38 


20 


-65 


25 


39 


36 


4.2 


2.97 


43.07 


^g2S 


16 


30 390.3 


Sy 1 


0, 


,0561 


18 


42 


09 


+79 


46 


17 


490" 


9.5 


2.0 


44.16 


< 0.1^'' 


23 


2E 1853.7+1534 


Sy 1 


0, 


,084 


18 


56 


00 


+ 15 


38 


13 


761 


9.3 


1.74 


44.48 


7 


10 


IH 1934-063 


Sy 1 


0, 


,010587 


19 


37 


33 


-06 


13 


05 


684 


2.7 


0.48 


42.08 


7 


4 


NGC 6814 


Sy 1.5 


0, 


,005214 


19 


42 


41 


-10 


19 


25 


488 


12.1 


2.92 


42.24 


< 0.05'' 


4 


IGR J19473+4452 


Sy 2 


0, 


,0539 


19 


47 


19 


+44 


49 


42 


969 


5.9 


1.05 


43.86 


11" 


11 


Cygnus A 


Sy 2 


0, 


,056075 


19 


59 


28 


+40 


44 


02 


1376 


21.6 


3.24 


44.39 


20^^ 


4 


MCG +04-48-002 


Sy 2 


0, 


,014206 


20 


28 


35 


+25 


44 


00 


187 


3.1 


1.10 


42.70 


7 


2 


4C +74.26 


AGN 


0, 


,104 


20 


42 


37 


+75 


08 


02 


72° 


1.9 


0.93 


42.35 


0.21-'' 


23 


MRK 509 


Sy 1 


0, 


,034397 


20 


44 


10 


-10 


43 


25 


73 


8.6 


4.66 


44.11 


< 0.01^ 


4 


GJF B21 16+818 


Sy 1 


0, 


.086 


21 


14 


01 


+82 


04 


48 


192" 


4.8 


1.8 


44.51 


< 0.1^ 


23 


IGR J21247+5058 


AGN 


0, 


.020 


21 


24 


39 


+50 


58 


26 


213 


11.9 


4.15 


43.57 


7 


20 


NGC 7172 


Sy 2 


0, 


.00868 


22 


02 


02 


-31 


52 


11 


401" 


15.1 


3.3 


47.74 


9.0^® 


23 


3C 454.3 


Blazar 


0, 


.859 


22 


53 


58 


+ 16 


08 


54 


92" 


2.8 


3.1 


46.56 


7 


26 


MR 2251-178 


Sy 1 


0, 


.063980 


22 


54 


06 


-17 


34 


55 


489 


7.0 


1.20 


44.07 


< 0.19^ 


4 


MCG -02-58-022 


Sy 1.5 


0, 


.046860 


23 


04 


44 


-08 


41 


09 


489 


3.9 


1.20 


43.79 


< 0.08^ 


4 



''IBIS/ISGRI exposure time 
'■flux /(2o-40keV) 10~ erg cm~^ s~ ^ 
^intrinsic absorption in 10^^ cm~^ 
^tentative redshift 

"not covered by survey presented here 

'references. — (1) Donato, Sambruna, Gliozzi 2005; (2) Bassani ct al. 2006; (3) Tartarus database; (4) Bcclimann ot al. 2006a: 
(5) Matt et al. 1997; (6) Lutz et al. 2004; (7) Bird et al. 2006; (8) Virani et al. 2005; (9) Plan ct al. 2005; (10) Masetti ct al. 2005; (11) 
Sazonov et al. 2005; (12) Beckmann et al. 2005; (13) Beckmann et al. 2004; (14) Courvoisier et al. 2003b; (15) Revnivtscv ct al. 2006; 
(16) Soldi et al. 2005; (17) Sazonov & Revnivtsev 2004; (18) Levenson, Weaver, & Heckman 2001; (19) Matsumoto, Nava, Maddox et 
al. 2004; (20) Masetti et al. 2004; (21) De Rosa et al. 2005; (22) Young et al. 2002; (23) This work; (24) Sazonov et al. 2004; (25) 
Revnivtsev et al. 2004; (26) Plan et al. 2006; (27) Gliozzi, Sambruna, Eracleous 2003; (28) Akylas, Georgantopoulos, Comastri 2001; 
(29) Piconcelli et al. 2006; (30) Kuiper et al. 2005 
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